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Abstract 

Non-perturbative parameters to the order 0( — ) in HQET expansion, A, Ai and A2, 
as well as the masses, m& and m c , are estimated phenomenologically from B and D 
meson system spectroscopy. We found relatively large value of nib and quite small 
value of A compared to previous estimates. 
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[I] The heavy quark effective theory (HQET) provides a systematic expansion in — 
to study physics of heavy hadrons containing a heavy quark [[[]], and has been applied to 
various area of phenomenology 0, [|. In mq — * oo limit the effective theory has a spin-flavor 
symmetry (heavy quark symmetry). Applied to the heavy hadron spectra, the symmetry 
leads to sum rules relating masses of heavy hadrons with different heavy flavors. E.g. using 
the sum rules, the masses of some beauty hadrons have been predicted in terms of the masses 
of known charm hadrons M . 

The observed mass splitting of the ground-state spin doublet is 141 (or 46) MeV for 
D (or B) mesons. It implies that the ^ corrections to the leading order are important 
for the realistic heavy hadron spectra, especially, for charm hadrons. Furthermore, the 
parameters of the order O(^) represent the contributions from the kinetic and chromo- 
magnetic operator || and are the fundamental parameters of HQET, which play a crucial 
role in the description of power corrections to the leading order. For example, the leading 
power corrections to B — > X q + ev decay rate are completely determined in terms of the 
parameters 0, which parameterize the nonperturbative effects. However, unfortunately, 
the parameters have not been calculated theoretically in an unambiguous way. Therefore, 
to estimate them from phenomenological analyses is quite meaningful. In this short note 
we study the spectra of heavy mesons to the order O(-) in the framework of heavy quark 
expansion, and estimate the parameters of the order 0{-^) from phenomenological analyses. 

From the effective Lagrangian in HQET, the mass of a heavy hadron can be written as 

M h = m Q + A + — (a+ < s Q ■ j t > b) , (1) 

where A is the "mass" of light degree of freedom from the binding energy and masses of the 
light quarks, a and b are the parameters which characterize the effects of the heavy flavor 
symmetry breaking and spin symmetry breaking at the order 0( — ), respectively: 

A = lim (M h -m Q ), 

mQ— »oo 

a = Al, 

2 ' 

b = 2A 2 . 

Evidently A, a and b do not depend on heavy quark mass, mq. In Eq. (1) sq is the heavy 
quark spin, and ji is the spin of the light system. The measurements of D and B meson 
masses are now quite accurate, and from Eq. (1) we have 

M D = -(M D + 3M D *) =m c + A + — ~ (1973 ± 2) MeV, (2) 
4 m r 



AM D = M D * - M D = — ~ (141 ± 1) MeV, (3) 

m c 

M B = ~(M B + 3M B *) =m b + A + — ~ (5313 ± 2) MeV, (4) 
4 m b 

A Mb = M B * - M B = — ^46 MeV, (5) 

"'/, 

where we also show the experimental values from PDG96 ||. 

Eqs. (2-5) have five unknown independent quantities, (m c , m b , A, a and b), with four 
independent observables. We choose m b as a free parameter^, so that we find the other 
parameters as functions of m b : 

m c = r x m b ~ 0.326 x m b , (6) 

A = (lM D + ~ +r ) x m b ~ ( 6928 ± 2 ) MeV ~ L326 x m ^ ( 7 ) 

a = rxm b -(m b - ^ 0.326 x m b ■ (m b - (4955 ± 3) MeV) , (8) 

b = AM B xm b ~ 46 MeV x m 6 , (9) 

where 

r = — = — — ^ ~ 0.326, 
m b AM D 

AM = M B -M D ^ (3340 ± 2) MeV. 

It is obvious that Eqs. (6-9) are not the definitions of the parameters, but phenomenological 
relations satisfying the experimental results of Eqs. (2-4). 

If we fix m;,=4800 MeV, which is deduced from a QCD analysis of the T system, we get 

m c ~ 1565 MeV, A ~ 563 MeV, a ~ -155 MeV x m c , b ~ 0.221 GeV 2 . 

Note that the value of a becomes negative. This nagaive value for a, certainly, can not be 
acceptable because the parameter a is essentially the expectation value of the kinetic energy 
operator of the heavy quark in the rest frame of the hadron in HQET H 0, and should be 
positive. This implies that the effective mass of the b quark inside the B meson is different, 
probably within order O(Aqcd), from that inside the T system if the expansions of heavy 
hadron masses in — are valid. 

In Fig. 1 we show the values of m c (in GeV), A (in GeV), a = — |Ai (in GeV 2 ) and 
b = 2X 2 (in GeV 2 ) as functions of m b . As shown in Eq. (8), only for m b > 4958 MeV the 

3 Here the parameter m b is phenomenologically defined only through Eq. (1) as an expansion parameter. 
It might be a pole mass or a running mass or else depending on given extra assumptions. In any case, 
however, it together with other parameters must satisfy Eq. (1) and the experimental constraints of Eqs. 
(2-5). 
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value of a becomes positive. Similary, only for m b < 5223 MeV the value of A can be positive 
from Eq. (7). The shaded region, 4958 < m b < 5223 (in MeV), represents for both a > 
and A > 0. From the shaded region we can find the region of parameter space allowed in 
HQET expansion to the order O(-): 

m b = [4958, 5223] MeV, (10) 
m c = [1616, 1703] MeV (or (m b -m c ) = [3342, 3520] MeV), (11) 



A = [0, 356] MeV, (12) 
a = --Ai = [0, 0.461] GeV 2 , (13) 
b = 2A 2 = [0.228, 0.240] GeV 2 . (14) 



[II] Previously the values of parameters have been estimated, e.g. see Ref. M, as: 

m b = (4.71 ±0.07) GeV, 
m b -m c = (3.39 ±0.04) GeV, 

A = (0.57 ±0.07) GeV, (see [8]) 
-Ai = (0.3 ±0.2) GeV 2 , (see [9, 10]) 
A 2 « 0.12 GeV 2 . 

Compared to our numerical bounds, Eqs. (10-14), previos estimations give relatively small 
value of m b and quite large A, however at the same time, very compatible values for (m b —m c ), 
Ai and A 2 . As is well known, the heavy quark mass difference, (m b — m c ), and spin symmetry 
breaking, A 2 , are precisely determined within 5 % from B and D meson spectroscopy. There 
are intrinsic uncertainties for the definitions of rriQ, A and Ai related to higher order pertur- 
bative corrections 0, 11]. However, we note that Eqs. (1-9) are perfectly valid independent 



of the definition of parameters in heavy quark mass expansion up to O(^). Therefore, the 
bounds on the parameters, Eqs. (10-14), and the experimental constraints, Eqs. (2-5), have 
to be satisfied simultaneously as pre-requirements: e.g. a parameter set chosen from central 
values of the previous estimates 

(m b , m c , A, - Ai, A 2 ) = (4.71 GeV, 1.32 GeV, 0.57 GeV, 0.3 GeV 2 , 0.12 GeV 2 ) 

can not simulataneously satisfy the experimental results, Eqs. (2-5), and should be discarded 
as a heavy quark expansion parameter set, even though an individual parameter can have 
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value outside or inside the bound depending on its proper definition. We also note that if 
the value of mj is chosen as ~ 4.7 GeV, the only parameter set consistent with Eqs. (2-5) is 
given as 

(m b , m c , A, - Ai, A 2 ) ~ (4.7 GeV, 1.5 GeV, 0.7 GeV, -0.8 GeV 2 , 0.11 GeV 2 ), 

which does not satisfy the bounds at all. 

Recently, it has been an important subject to obtain an accurate value of the kinetic 
energy, /z 2 (= — Ai), of the heavy quark inside 5-meson. Ball et al. |3| calculated using 
the QCD sum rule approach and obtained /i 2 ~ 0.50 GeV 2 for 5-meson, while Neubert 
M obtained — Ai ~ 0.1 GeV 2 . It should be noted that those two derivations differ in the 
choice of the 3-point correlation functions used to estimate the matrix elements of interest. 
Bigi et al. derived an inequality between the expectation value of the kinetic energy 
operator of the heavy quark inside the hadron and that of the chromomagnetic operator, 
A* 2 > |(My 2 — Mp 2 ), which gives /i 2 > 0.36 GeV 2 for S-meson system. However, Kapustin 
et al. || showed later that this lower bound could be significantly weakened by higher order 
perturbative corrections to the 3-point functions. Hwang et al. || also calculated the value 
by applying the variational method to the relativistic Hamiltonian, and obtained /U 2 ~ 0.44 
GeV 2 . Similarly de Fazio || computed the matrix elements of the kinetic energy operator 
by means of a QCD relativistic potential model, and found // 2 ~ 0.46 GeV 2 . Besides the 
theoretical calculations of /i 2 , Gremm et al. [[TI| extracted the average kinetic energy by 
comparing the prediction of the HQET [[| with the shape of the inclusive B — > Xlv lepton 
energy spectrum, and obtained —\\ = 0.19 ± 0.10 GeV 2 . 

If we choose the previously estimated average Ai value [0, £J [TIJ, — Ai = (0.3 ±0.2) GeV 2 , 
we can find quite narrow region of parameter space from the given experimental bounds, 
Eqs. (2-9), and Fig. 1: 

m b = (5.05 ± 0.07) GeV, (m b - m c ) = (3.41 ± 0.04) GeV, A = (0.23 ± 0.08) GeV. (15) 

1 /5 

The experimental value of semileptonic decay width of B meson can give the value rn b riQ CD 
12j| , with QCD correction t/qcd, by comparing with theoretical inclusive semileptonic decay 



width. After using B s \ and t b from PDG96 and using t/qcd = 0.77 ± 0.05 [|13|], we can 
get m b = (5.1 ± 0.2) GeV, which is in good agreement with our result. It is interesting to 
note that the value of hard pole mass m b obtained by Chernyak [TTJI is very close to that in 



Eq. (15). Finally we note that higher order corrections of 0(-^r) can change the bounds of 



Q 



our estimation. However, as shown in [|, the correction to m b is about ~ ±0.004 GeV, 
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and we can safely conclude that the higher order effect can not affect our results at all. 
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Figure 1: The values of m c (in GcV), A (in GeV), a = (in GcV 2 ) and b = 2A 2 (in GeV 2 ) as functions 

of m b . The shaded region, 4958 < m b < 5223 (in MeV), represents for both A > (sec Eq. (7)) and a > 
(see Eq. (8)). 
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